CONTRIBUTIONS TO THE THEORY OF THE 
BARNES FUNCTION 

V. S. ADAMCHIK 

Abstract. This paper presents a family of new integral represen- 
tations and asymptotic series of the multiple gamma function. The 
numerical schemes for high-precision computation of the Barnes 
gamma function and Glaisher's constant are also discussed. 



1. Introduction 

In a sequence of papers published between 1899-1904, Barnes intro- 
duced and studied (see j3 CHI EEU 112] ) a generalization of the classi- 
cal Euler gamma function, called the multiple gamma function T n (z). 
The function T n (z) satisfies the following recurrence-functional equa- 
tion HHJ|32|: 



(1) 



r n+1 (^ + i) = ^±i^, zee, nen, 

Uz) = T(z), 

r„(i) = i, 



where T(z) is the Euler gamma function. The system (0) has a unique 
solution if we add the condition of convexity 



Jn+1 

(-ir +1 ^iogr n (z)>o, z>o. 
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In this paper we aim at the most interesting case G(z) = 1/T 2 (z), 
being the so-called double gamma function or the Barnes G-function: 

G{z + l)=T{z)G{z), zeC, 
G(l) = G{2) = G(3) = 1, 
d 3 

— \ogG(z)>0, z>0. 

dz A 

The G-function has several equivalent forms including the Weier- 
strass canonical product 

G( , + 1) . w i exp( _£±£!^))n( 1 + £)^ P (i!-,) 

k=l 

and the functional relationship with the Hurwitz zeta function ((t, z), 
due to Gosper 20J and Vardi [31J 



(2) logG(z + l)-^lo g r(z) = C'(-l)-C / (-l^), &(z)>0, 

where 7 denotes the Euler-Mascheroni constant and the derivative 
('(t, z) = 4C(*> z ) is defined as an analytic continuation provided by 
the Hermite integral (j!8|) : 

9 2 
z z z 

C'(-M) = y log ^ - — - -logz + 

(3) 

f°° arctan(^) , f°° x log(x 2 + z 2 ) , 







e""'-" - 1 ./n e 



In particular, 

(4) f <(_ 1) = <<(_!, ) = 2 jTji^dx. 

For z — > 00, the G-function has the Stirling asymptotic expansion 

log G(z + 1) = I log(27r) - log A - ^- + ^ - -L) log * + O (I 
where A is the Glaisher-Kinkelin constant given by 

(5) logA = l-C'(-l). 
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Combining this definition with the representation (J3J), we obtain a new 
integral for log A: 

(6) hw i r^ fc 

v ; 6 12 2tt 2 J e x - 1 

The constant A originally appeared in papers by Kinkelin [23] arid 
Glaisher ^7J El] on the asymptotic expansion (when n — > oo) of 
the following product 

1 1P 2 2P ...n" P , peN. 

Similar to the Euler gamma function, the G-function satisfies the mul- 
tiplication formula |31| : 

G(n Z) = e C'(-l)(l-« 2 V 2 * 2 /2-n, + 5/12 

flj\ 71-1 71-1 . . 

[) (2^)^ n2) nn G (^+^)' nen 

The Barnes function is not listed in the tables of the most well-known 
special functions. However, it is cited in the exercises by Whittaker and 
Watson 33, and in entries 6.441, 8.333 by Gradshteyn and Ryzhik 
|2Tj . Recently, Vardi [3*T] . Quine and Choi j2ZJ, and Kumagai [21] 
computed the functional determinants of Laplacians of the n-sphere 
in terms of the r n -function. Choi, Srivastava [T3*l ITi] and Adamchik 
j3 HE] presented relationships (including integrals and series) between 
other mathematical functions and the multiple gamma function. 

The T n function has several interesting applications in pure and ap- 
plied mathematics and theoretical physics. The most intriguing one 
is the application of G-function to the Riemann Hypothesis. Mont- 
gomery j2H] and Sarnak [21] (see [22] for additional references) have 
conjectured that the non-trivial zeros of the Riemann zeta function are 
pairwise distributed like eigenvalues of matrices in the Gaussian Uni- 
tary Ensemble (GUE) of random matrix theory. It has been shown in 
works by Mehta, Sarnak, Conrey, Keating, and Snaith that a closed 
representation for statistical averages over GUE of N x N unitary ma- 
trices, when N —>■ oo can be expressed in terms of the Barnes functions. 

Another interesting appearance of the G-function is in its relation to 
the determinants of the Hankel matrices. Consider a matrix M of the 
Bell numbers: 
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M„ 



( B\ B 2 ... B n \ 
B 2 B 3 ... B n+ i 



\ B n B n+ i ... B 2n -i J 
where B n are the Bell numbers, defined by 



k=0 

and {^} are the Stirling subset numbers 
n — ll [n — 1 



k 



k 



k-1 



1, n = 0, 
0, n ^ 



It was proved in j^J, that 

detM n = G{n + 1), nGN. 



2. Special cases of the Barnes Function 

For some particular arguments, the Barnes function can be expressed 
in terms of the known special functions and constants. The simplest 
(but not-trivial) special case is due to Barnes 0: 



(9) 



logG 



1 log2 logTT 3 

- H log A 

8 24 4 2 6 ' 



where A is the Glaisher-Kinkelin constant (JHJ). The next two special 
cases are due to Srivastava and Choi |14j : 

(io) logo(-) =1-9--* \o g r(-) - -log a, 



. ogG (5)=io gG (I) + | 



G log 2 log Ti 
2t7 8 4~ 
where G denotes Catalan's constant: 



iogr(i 



G-V ( " X 



k=0 



(2k + iy 



Identities (fTU|) and (|TT|) are proven by using the duplication formula 
((Zj) (setting n = 2) together with (jHJ). More special cases for z — |, ^, | 
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and | can be derived as well. Here is one of such formulas (see j3] for 
the others): 

Proposition 2.1. 



.37 9 72 18V3 

-logr - --log a ^ uy 



3*^ 3" b " 12^' 



where ijj^ l \z) = ^logT(z) is the poly gamma function. 



Proof. We recall the Lerch functional equation for the Hurwitz zeta 
function (see [21 IE]) 

C(-n, z) + (-l)T(-^ 1 - «) = — — r B„+i(z) + 

n + 1 

e-™ /2 ^Li n+1 ( e 2 ^), < * < 1, nGN, 

where z) = -^((t,z), B n (z) are the Bernoulli polynomials and 
Li n (z) is the polylogarithm. Setting n = 1 and using (J2J) yields the 
reflexion formula for the G-function (see also [TH Hb|): 

where Ci2(z) is the Clausen function defined by 



(14) Cl 2 (z) = -9f(Li 2 (e- fa )). 

Upon setting 2 = | to (fTB^l. and making use of 

cl2 (2^_V> (1) (!) 27T 2 



3 / 3^ 9^' 
we obtain 
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On the other hand, using the multiplication formula (JJJ) with 
and n = 3, we find 



z 



(16) GilU 2 " 37A2e2/3 



3 



.3/ V3/ y27rA 8 r(|)' 

Now, combining (fTKJ) and (jlfjjl . we conclude the proof. □ 

It remains to be seen if the Barnes function (or the multiple gamma 
function) can be expressed as a finite combination of elementary func- 
tions and constants for other rational arguments. The most general 
result available for the multiple gamma function T n (z) is a closed form 



at z — 2 • 



Proposition 2.2. For n E N 

(17) 



i)>-i)!io g r,(i) = - < 2 "-f°s% 



k=l v ; k=l 



where Pk, n {\) are coefficients by x k , k < n in expansion of 



n—l j 

H( x + 3- o) 



Proof. The proof follows directly from formula (17), by setting 
z = \. □ 

Here few particular cases: 

the gamma function (n = 1): 

r(i) = ^ 

the Barnes function (n = 2): 

1 2 l/24 e l/8 



'2 7 A 3 / 2 7rV4 

the triple gamma function (n = 3): 

r .l. A 3/2 ^ ^7C(3) 1 

r 3(o) = o1/ M eX P 



'2 7 2V24 ^ \ 32 7T 
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3. Hermite integrals 

In this section we derive a few integral representations for the Barnes 
function. We start with recalling the Hermite integral (see [111213) for 
the Hurwitz zeta function: 

, . . z~ s z 1 ^ 3 f°° sin(s arctan(a;/z)) 

(18) C(s.z) = + + 2/ i \J—Lt dx, 

which provides analytic continuation of ((s, z) to the domain s G C — 
{1}. Differentiating both sides of (fT8|) with respect to s, letting s = — 1 
(and s = —2), and using the second Binet formula for logT(z) (see |Zj) 

/ \ , , i . 1., log(27r) f°° arctan(x/z) , 

(19) \ og T(z) = (z--)\ogz-z + ^ r ^ + 2j^ —^-L-ldx 

together with (J2J), we readily obtain: 

2 

logG(z + 1) = y (log^ - H 2 ) - (zC(0) - C'(-l)) - 

(20) 

" x xlog(x 2 + z 2 ^ 



e 2irx _ I 



dx, $l(z) > 0, 



where H p = Xwb=i are the harmonic numbers and ('{z) is a deriv- 
ative of the Riemann zeta function. 

In the same manner, repeating the above steps n times, we derive 
the general integral representations. 

Proposition 3.1. For $i,(z) > 
(2n)! logT 2n+1 (z + l) = 



z 2n+l(l 0gz _ H2n+l) - / 2n 



2n + 1 

k=0 



(21) 



E(-i) fe fc! ' io g r fc+1 (z + i) + 

fc=l ^ J 



oo ^2n 



arctan (x/ z) 



s 
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and 



(2 n -i)!io E r 2 „(, + i) = - ^ (log 2 ;-^ ) + 

E(- 1 )'( 2n r 1 )<'(-*> a "-*- 1 + 

! — n \ / 



(22) 



fc=0 
2rc-2 



fe=i ^ J 



e 27rs - 1 

where j^} are £/ie Stirling subset numbers, defined in (0). 

Upon substituting 2 = 1 into formulas (}2*Tj) and (j22J), we obtain new 
integrals 

(23) 2 (-D» [ + 1 (-D*( 2 ;)c<(-*) 

(24) 

17 u fc=0 v 7 

Performing simple transformations of the integrands in (J23J) and 
will apparently lead to a large number of special integrals: 



DC 







, x log X 

2 I ^rzi dx = ?(-i) 



/°°J^^=12C'(-D + log2 

r»log(l + ^) 3 23 1 
J e 2 ™ + l 4 24 6 2 s v ; 



r~>o 



, arctan (-) log(27r) , 1 N s 

2 I ^ dx = z\ogz-z + -*±-L - \ogT(z + -), R(z)>0 



'0 
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4jf ^gU=-2 + 31og2 



-= ?rr = log 2: — ib(z) , argfz) ^ — 

- l)(x 2 + z 2 ) & ^ y 1 2z SV ; T 2 



4. BlNET-LIKE REPRESENTATION 

In this section we derive the Binet integral representation for the 
Barnes function. 

Proposition 4.1. The Barnes G-function admits the Binet integral 
representation: 

2 i 

\ogG(z + 1) = z\ogT{z) + Z --^L b 2 {z) - log A - 

r°° e~ zx / 1 1 1 x \ 
J ar VI — e x x 2 12/ 

Proof. Recall the well-known integral for the Hurwitz function |H] 
(26) C(s, z) = TvT / 1 6 a ^ > X > ^) > 

r(s) y i - e * 

The integral in the right hand-side of (j2T?j) can be analytically continued 
to the larger domain 3?(s) > —2 by subtracting the truncated Taylor 
series of = 1 _ x at x — 0. Since 

and 



T(s) ./n Vx ' 2 ' 12/ ~ 12 " 2 1 - s' 

&(s) > > 

we therefore can rewrite (l26l) as 
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z 1 s Z o 



t(s,z) = 1 1 z~ s ~ y + 

sv ' ; 1 - s 2 12 



/°° x^V** f - - - — )dx, > -2. 

r(s) y Vl-e-* x 2 12 J ' v ; 

Differentiating the above formula with respect to s and computing the 
limit at s — —1, we obtain 



C'(-M) = --- + ^fl a (*)+ 

/•OO -ZX , 1 1 1 . 

/ —(—-"2-12)^' *w>°- 

where B 2 (z) is the second Bernoulli polynomial. The convergence of 
the above integral is ensured by (J27)l . Finally, using (J2J), we conclude 
the proof. □ 



5. Asymptotic expansion of Y n {z) 

The asymptotic expansion for T n (z) when z — > 00 follows straight- 
forwardly from formulas (|21|) and ()22|) by expanding log(l + x 2 ) and 
arctan(x) into the Taylor series, and then performing formal term-by- 
term integration. 

Lemma 5.1. For 9?(z) > 



log G(z + 1) = y (log z-^-^-z C(0) + C'(-l) 

( 28 ) 

E B2k+2 n( 1 x 



Proof. We start with the integral representation (|20p. By recalling the 
Taylor series for log(l + x) 

X ™ C_ \)k x k+1 (—l) n f x y n+1 

Mi + -) = E tzt (-) -^/ V dy. 

z ^ k+1 z z n+1 J y + z 

we find that 
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x\og(l + X 2 Z 2 ) . 
o dx 

e 27TX _ I 



i r°° x 2k+3 _ (-i) n r°° tdt r t2 y n+1 

/L k + l z 2k+2 J o e 2nx _ I X ~ ~^+2 J o e 2nt _ \ J q ^^2 



which, by appealing to (formula 1.6.4 in 

roc ±2k-l r> 

(29) / dt = (-l) fc+1 ^ 

reduces to 



x\og(l + x 2 /z 2 ) 
n ; dx 

e 2nx _ I 



y B2 k+ 2 . )n r tdt [*/' y W 

f-^4k(k + l)z 2k 1 ' J e M -l i 1+y 2 



Substituting this into (|2(jp . we complete the proof. 



Lemma 5.2. For 9R(z) > 

2 log T 3 (z + 1) = -j ( log z - ^) + y ( log z - | + C'(0)) - 

^ - z(C(0) + 2C'(-1)) + C'(-l) + C'(-2) - 
(30) „ 

E-£>2fc+2 \ ^ -D2fc+2 

4fc(fc + l)z 2fc ~ ^ 2 (jfe + 1) (2A; - 1) z 2 ^ 1 
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Proof. Employing the same technique as in getting (|28jl , we begin with 
the integral representation of the triple gamma function: 



2' 3 



(31) 



21o g r 3 (z + 1) = -^-(logz - H 3 ) + \ogG(z + 1) + 



^ 2 C / (0)-2<(-l) + C , (-2))+2 



00 x 2 arctan(x/ z) 



o 



dx 



that follows directly from Proposition 3.1 with n — 1. Expanding 
arctan into the Taylor series 

arctan(^) = £ ^-j-(-) + ^ / d V 
and making use of (|29|). we readily obtain 

/■°° x 2 arctan(x/z) _ i?2fc+2 
7 e 2 ™ - 1 = ~ ^ 4(fc+ l)(2fc- l)z 2fc -! 4 



fc= -^4(fc + l)(2A;-l)z 2fc - 1 
-l) n r 00 t 2 dt [* y 2n 



~2n— 1 / e 2 7rt 1 / ,,2 I ~2 

z Jo e — L Jo y + z 

The conclusion of the lemma follows immediately by substituting this 
into formula ()31j) and using Lemma 5.1. 

□ 



6. Implementation 

In this section we discuss a numeric computational scheme for the 
Barnes function. Let us consider an integral representation dated back 
to Alexejewsky jS] and Barnes [Oj: 

/ logT{x)dx = Z{1 ~ ^ + Z - log(27r) + zlogr(^) - logG(^ + 1), 
Jo 11 

where dt(z) > — 1. Employing integration by parts, we derive 

(32) log G{z+1) = ^—A + - log(27r) + f x ip{x) dx, > -1 

2 2 J 



The integral representation (j52|) provides an efficient numeric procedure 
for arbitrary precision evaluation of the G-function. This representa- 
tion demonstrates that the complexity of computing G(z) depends at 
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most on the computational complexity of the polygamma function. 
The later can be numerically computed by the Spouge approximation 
p?Uj (a modification of the Lanczos approximation). For the numerical 
integration in ()32|) we shall use the Gaussian quadrature scheme. 

The restriction > — 1 in (J32j) can be easily removed by analyt- 
icity of the polygamma. For example, by resolving the singularity of 
the integrand at the pole x = —1, we continue logG(z + 1) to the wider 
area dt(z) > -2: 

logG(z+l) = Z{ " l ~ Z) + Z - log(27r)+log(2;+l)+^ (x^(x) - dx 

Another method of the analytic continuation of the G-function is to 
use the following well-known identity 

ip(x) = if)(—x) — 7rcot(7rx) 

X 

which, upon substituting it into J32j), yields 



(33) logG(l-2) =logG(l + ^)-2log(2vr)+ / vrx cot(vrx) dx. 

Jo 

The identity holds everywhere in a complex plane of z, except the 
real axes, where the integrand has simple poles. Therefore, combining 
(|3*2*|) with (}3*3*|) yields the following definition of the Barnes function 

G(-n) = 0, nGN, 

(34) G(z) = (2tt)( z - 1 )/ 2 exp (- ^ - ^ - 2) + J x ^(x) dx 

arg(z) ^ 7T, 

where the contour of integration 7 is a line between and z—\ that does 
not cross the negative real axis; for example, 7 could be the following 
path {0, i, i + z — 1, z — 1}. 
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It remains to define G(z) for z e R~. We do this by using the 
reflexion formula (|13j) . which upon periodicity of the Clausen function 
(j!4j) can be rewritten as follows: 

G (-z) = (-1)^-^ + 2) Z+ \ 

(35) 

exp(i-Cl 2 (27r(^- L^J))), ZER-. 

7. Glaisher's constant 

In this section we discuss numeric computational schemes for Glaisher's 
constant: 

(36) A = exp(^-C / (-l)) 

In light of the functional equation for the zeta function 

CM = 2-7^^(1 - a )C(l - s) sin (^) 
we can rewrite (|3T)j) as 

(37) ^esp^-^)^. 

Exactly this formula is used in Mathematical '4. 2 for computing Glaisher's 
constant. 

Another computational scheme for Glaisher's constant is based on 
the Barnes G-function identity: 

(38) logA=l + ^-^-'logG(iY 
K ' 5 12 36 6 3 S \2J 

Here can be computed via integral (J3"2j) . 

A more efficient algorithm is provided by the following identity 

k=l 

which allows to approximate Glaisher's constant as 
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k=l 



This approximation is suitable for arbitrary precision computation - it 
requires |"| log 2 10] terms to achieve p decimal-digit accuracy. 

We close this section with a proof of formula Recall that jH] 



(40) COO - 1 = F7T / w z ^ > 1. 



T(s) Jo e rr (e 1 ' - 1) 

Upon replacing ((2k + 1) — 1 in by its integral representation (|3Uj) 
and inverting the order of summation and integration, we obtain 



. , , , log2 1 f 00 dx ^ t lk ( 3 6 
(41)logA = — +~l - I _ Ty g_(28 + — - — 

Next, noting that 

\ (2^) 



00 ^2fc 



and 



g (2k)\(k + s) = &] Q x2S '^ cosh ^ ' !) ^ 
we compute the inner sum 

j.4 00 +2fe / o c . 

L y * f 28 + -i- - -M = -36 + 3 1 2 - 14 t 4 + 
2^(2fc)!V 1 + k 2 + kJ 

(36 + 15 1 2 + 14 1 4 ) cosh(t) - (36 1 + 3t 3 ) sinh(t) 

In the next step, we substitute this back to (j4~Tj) and, employed by (J2HJ), 
evaluate the integral. Unfortunately, integration cannot be done term 
by term. The way around is to multiply the integrand by t x , A > 
and consider the limit as A — > 0. This yields 
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36t A ~ 4 36 t A " 4 72 t x ~ 4 36 t x ~ 3 

e* - 1 + e 2t (e* - 1) ~ e* (e* - 1) ~ e* - 1 + 

36 t A " 3 15t A ~ 2 15t A ~ 2 6t A ~ 2 3t A - x 

+ ^ ^ + ^TT^ 7T + ^^ 7T-- 7 + 



e 2 *(e 4 -l) e*-l e 2i (e* - 1) e* (e* - 1) e* - 1 
3t A-1 14 t A 14 t A 28 t A 

+ ~ + ~~KZ — "i — 7 7T - — 7T dt 



e 2 *(e*-l) e'-l e 2 * (e* - 1) e* (e* - 1) 
36T(A- 3) - 9 2 5 - A T(A-3) -36T(A- 2) - 9 2 4 - A T(A-2) - 

21 T(A - 1) - 15 2^ A T(A - 1) - 3 T(A) - + 14 T(A + 1) - 

ZE^-ti)+36r(A-l)C(A-l). 

Computing the limit as A — > 0, we find that the right hand-side of the 
above expression evaluates to 

3-31og2-36C'(-l). 
This completes the proof of (j39|) . 
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